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The first stars and quasars are known sources of liard ionizing radiation in the 
first billion years of the Universe. We examine the joint effects of X-rays and hard 
UV radiation from such first-light sources on the hydrogen and helium reionization 
of the intergalactic medium (IGM) at early times, and the associated heating. We 
study the growth and evolution of individual H II, He II and He III regions around 
early galaxies with first stars and/or QSO populations. We find that in the presence 
of helium-ionizing radiation. X-rays may not dominate the ionization and thermal 
history of the IGM at z ^ 10-20, contributing relatively modest increases to IGM 
ionization, and heating up to ^ 10^-10^ K in IGM temperatures. We also calculate 
the 21 cm signal expected from a number of scenarios with metal-free starbursts and 
quasars in varying combinations and masses at these redshifts. The peak values for 
the spin temperature reach ^ 10^-10^ K in such cases. The maximum values for the 
21 cm brightness temperature are around 30-40 mK in emission, while the net values 
of the 21 cm absorption signal range from ^ a few to 60 mK on scales of 0.01-1 Mpc. 
We find that the 21 cm signature of X-ray versus UV ionization could be distinct, with 
the emission signal expected from X-rays alone occurring at smaller scales than that 
from UV radiation, resulting from the inherently different spatial scales at which X-ray 
and UV ionization/heating manifest. This difference is time-dependent, and becomes 
harder to distinguish with an increasing X-ray contribution to the total ionizing photon 
production. Such differing scale-dependent contributions from X-ray and UV photons 
may therefore "blur" the 21 cm signature of the percolation of ionized bubbles around 
early halos (depending on whether a cosmic X-ray or UV background built up first), 
and affect the interpretation of 21 cm data constraints on reionization. 

Key words: stars: Population HI. galaxies: high-redshift. (galaxies:) quasars: general, 
galaxies: star formation, cosmology: theory, (cosmology:) dark ages, reionization, first 
stars. 



1 INTRODUCTION 

The first billion years after the Big Bang represents a period 
of great interest for studies of both galaxy formation and 
the evolution of the Universe as a whole. This period sees 
the formation of the first galaxies (Wise et al. 2008) and, 
consequently, the beginning and completion of the process 
of reionizing the Universe (Loeb & Barkana 2001; Loeb 2009) 
as a result of the copious number of ionizing photons emitted 
by these sources. Current and future facilities aim to probe 
this epoch of the Universe both using traditional methods 
such as surveying faint galaxies (e.g. the James Webb Space 
Telescope; Gardner et al. 2009) and using novel techniques 
such as 21cm cosmology (Furlanetto et al. 2006) to probe 



the distribution of neutral hydrogen during the process of 
reionization. Understanding this epoch of the Universe from 
a theoretical perspective therefore requires an understanding 
both of the sources of ionizing photons and of the thermal 
and ionization state of the intergalactic medium (IGM) at 
these times. 

Additionally, the thermal and ionization history of the 
IGM as a function of cosmic redshift, z, strongly aflects the 
"visibility" of the most distant galaxies and quasars (Madau 
1995; Meiksin 2006; Dayal et al. 2011), and the feedback ex- 
erted on the formation of new galaxies (Efstathiou 1992; 
Quinn et al. 1996; Navarro & Steinmetz 1997; Barkana & 
Loeb 1999; Bullock et al. 2000; Somerville 2002; Benson 
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et al. 2002a,b; Koposov et al. 2009; Muoz et al. 2009; Busha 
et al. 2010; Macci et al. 2010). The process of reionization 
is expected to begin with the formation of ionized bubbles 
around luminous sources in the redshift range z — 10-20. 
These bubbles will eventually grow in size and number until 
complete overlap is reached and the Universe becomes fully 
reionized. The shapes and sizes of bubbles will be controlled 
by the cosmological density field and the process of galaxy 
formation. Their internal ionization and temperature struc- 
ture will depend on the spectrum of the input source (i.e. 
how hard the photons are) and the efficiencies of recombi- 
nation and cooling processes. 

The recent data from WMAP-7 (Larson et al. 2011) re- 
veal that the IGM is fully ionized up to 2; ~ 10, most likely 
with a period of partial ionization at higher redshifts. Theo- 
retical work over the last fifteen years has focused mostly on 
the hydrogen reionization of the IGM (Gnedin &: Ostriker 
1997; Chin & Ostriker 2000; Ciardi et al. 2000; Somerville 
et al. 2003; Onken & Miralda-Escude 2004; Benson et al. 
2006; Furlanetto et al. 2006). However, helium reioniza- 
tion has received comparatively less attention, ranging from 
calculations of helium/hydrogen reionization from the first 
stars and QSOs at z > 6 (Venkatesan et al. 2003; Wyithe 
& Loeb 2003) to studies of helium reionization by QSOs 
at z ~ 3 (Sokasian et al. 2003; Furlanetto & Oh 2008a,b). 
Although helium is the second most abundant element, its 
substantially higher ionization energy relative to hydrogen, 
as well as its interactions with X-rays through secondary ion- 
izations, can lead to significant effects for the high-z IGM 
and the cosmic microwave background (CMB) once reion- 
ization has occurred even to a partial degree. Additionally, 
X-rays have greater penetrating power relative to UV radia- 
tion. When occurring in combination with helium ionization 
from the first stars and quasars. X-rays could act to strongly 
alter the ionization and thermal history of the IGM. 

In this work we investigate the joint impact of X-rays 
and helium-ionizing radiation from the first galaxies on IGM 
reionization and heating. We focus on the growth and evolu- 
tion of individual ionization fronts in H and He, rather than 
a fully evolving cosmological calculation, which we plan to 
pursue in future work (§4). We study whether the differ- 
ing contributions arising from X-rays versus UV ionization 
can be distinguished through 21 cm observations. Recent 
papers by other authors have focused on specific aspects 
of this problem in other contexts, e.g., helium reionization 
by quasars at lower redshifts (z ~ 3) (Bolton et al. 2009; 
McQuinn et al. 2009; Furlanetto & Oh 2008b), without ex- 
plicitly considering the effects of X-ray heating (Furlanetto 
& Oh 2008a), or, with only a single high-mass star embed- 
ded in a high-z galaxy halo (Chen & Miralda-Escude 2008). 
We will demonstrate that X-rays may not play a dominant 
role in high-z ionization, contrary to the theoretical expec- 
tations in some previous works (see, e.g., Thomas & Zaroubi 
2008), and will ask the question: does there exist a cosmolog- 
ical epoch when the IGM's thermodynamic and ionization 
properties are determined mostly by X-rays? 

The tradeoff of these ionization effects will have impor- 
tant consequences for predictions for future radio observa- 
tions that plan to see ionized bubbles in emission or absorp- 
tion against the CMB. There has already been a substan- 
tial body of work on the feedback on ionization (Venkatesan 
et al. 2003; Tumlinson et al. 2003; Wyithe & Loeb 2003) and 



emission line signatures (Oh et al. 2001; Tumlinson et al. 
2001; Venkatesan et al. 2003; Dawson et al. 2007) arising 
from first-light sources that have hard ionizing spectra. Here, 
we focus on the radio signatures as the topology of reion- 
ization arising from X-rays versus UV radiation is expected 
to be different. We also test other theoretical predictions for 
the growth of individual ionized regions around early galax- 
ies, e.g., that for sufficiently hard sources such as the first 
stars and QSOs, the H and He I-fronts may track each other 
closely. 

The remainder of this paper is arranged as follows. In 
Section 2, we describe the model that we use to follow the 
growth of cosmological ionization fronts around evolving 
sources. In Section 3, we present our results for the thermal 
and ionization properties of such regions, their observable 
signatures (including 21 cm signals) around a set of rep- 
resentative sources, and compare our findings with earlier 
works in this field. We conclude in Section 4. 



2 BACKGROUND AND MODELS 

We assume a background cosmology using the most recent 
cosmological parameters fits from the WMAP-7 CMB data 
(Larson et al. 2011). We combine the formalism for study- 
ing the non-equilibrium evolution of hydrogen and helium 
in the IGM in Venkatesan et al. (2003) and Tumlinson et al. 
(2004), and the input processes related to X-ray ionization 
in the high-z IGM in Venkatesan et al. (2001) with the code 
Galacticus. Galacticus is a newly developed semi-analytic 
code on galaxy formation (Benson 2011) which includes feed- 
back from high-redshift star/quasar formation while meeting 
current experimental constraints at lower redshifts. Here, we 
have utilized it to solve for the growth of a spherical ioniza- 
tion front around a point source in the IGM. The ionizing 
and heating processes included in this code are described in 
detail below. 

We are primarily interested in the effects of hard ioniz- 
ing radiation from the first galaxies - these are assumed to 
be of order 10* Mq in total mass and of approximate size 
1 to a few kpc^. We follow the advancing ionized fractions^ 
around a starburst and/or quasar in such a halo, and treat 
the IGM as being homogenous around the source. In partic- 
ular, we do not include a density enhancement as would be 
expected if the source forms in the center of a dark matter 
halo. In general, the ionization fronts we find are much larger 
than the sizes of typical halos at these redshifts and so will 
be insensitive to the details of the density profile on small 
scales. Additionally, sources such as those considered here 
will likely form in halos sufficiently massive to coUisionally 
ionize hydrogen and helium, such that the photoionization 



^ A lO'^'^ M0 halo at z = 10 has an approximate physical (not 
conioving) virial radius of ^ 7 kpc, with a galaxy of size ^ 1 kpc 
in it. A massive Milky Way-mass halo (10^^ M0) at that redshift 
would be about 33 kpc, with a typical galaxy of a few kpc in size. 
■^ This is in contrast to Venkatesan et al. (2001), where the av- 
erage IGM ionization fraction was computed without tracking 
the growth of individual I-fronts around the halos containing the 
QSO. 
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front would begin growing from the edge of the collisionally 
ionized region^. 

We consider quasars with varying black hole (BH) 
masses, and model a typical QSO spectrum with the fit given 
in Haardt & Madau (1996). We assume that the duty cy- 
cle of the QSO is 100 million years typically - significantly 
longer duty cycles would exceed the Hubble time at z ~ 
10-20. In our models, we allow the AGN to be on for 100 
Myr before it is shut off. We include the efi'ects of metal-free 
stars occurring in starbursts of varying masses — the fits 
are taken from Venkatesan et al. (2003). 

The non-equilibrium ionization fractions are calculated 
including the following processes: photoionization, coUi- 
sional ionization, case B radiative recombination, dielec- 
tronic recombination for He I, and the coupling between H 
and He caused by the radiation fields from the He I 24.6 eV 
recombination continuum and from the bound-bound tran- 
sitions of He I (Venkatesan et al. 2001). The photoionization 
cross sections for H I and and He II are taken from Spitzer 
(1978), and from Verner et al. (1996) for He I. The ratio of 
the H I to He I photoionization cross sections decreases with 
photon energy, ranging from about 5% at 100 eV to 3.5% at 
1 keV. This implies that an X-ray photon is "seen" better 
by a He I atom than by a H I atom. 

We also include secondary ionizations and excitations 
of H I and He I arising from the X-rays (ShuU & van Steen- 
berg 1985). As noted in Venkatesan et al. (2001), a typical 
X-ray photon is far more likely to be absorbed by He I rather 
than H I, so that secondary ionization (rather than direct 
photoionization) is most relevant for H I when X-rays domi- 
nate photoionization. The resulting photoelectrons will ion- 
ize many more H I atoms than He I, H I atoms being more 
numerous. As the background ionization increases, the pho- 
toelectron deposits more and more of its energy in heat and 
less in collisional ionizations/excitations. ShuU & van Steen- 
berg (1985) assumed that the ionization fractions of H I and 
He I were equal, and we have replaced the generic ioniza- 
tion fraction in their formulae with the electron fraction x^, 
which is more directly relevant for the IGM. 

The thermal evolution of the gas is computed including 
the following processes (Venkatesan et al. 2001): photoelec- 
tric heating from the secondary electrons of H and He, which 
is itself a function of the background ionization levels (Shull 
& van Steenberg 1985), and, heating from the H I photo- 
electrons liberated by the bound-bound transitions or the 
24.6 eV recombination continuum of He I. Cooling terms in- 
clude radiative and dielectronic recombination (Venkatesan 
et al. 2001 and references therein), thermal bremsstrahlung, 
Compton scattering off the CMB, collisional ionization and 
excitation, and the adiabatic expansion of the IGM. The 
contributions to heating and cooling from the scattering of 
the secondary Lya photons from X-ray ionization is negligi- 



^ In a fully 3-D calculation these halos would accrete most of their 
mass via cold filaments of gas which are not shock heated as they 
enter the halo and so are not collisionally ionized. It is beyond 
the scope of this work to examine the effect of such filaments on 
the growth of ionization fronts (Keres et al. 2005), but they can 
be expected to impede the growth of the front along directions 
coinciding with a filament, while permitting faster growth along 
directions between filaments. 



ble (Chen & Miralda-Escude 2004, 2008) and is not included 
here. 

Our ID non-equilibrium ionization code includes all of 
the above ionization and heating processes, and solves for 
the evolution of the thermal and ionization state around the 
source as follows. The IGM surrounding the source is divided 
up into a large number of concentric spherical shells. Unless 
otherwise noted, we use 1000 shells, spaced logarithmically 
in radius from 10~^ to 10 Mpc. These shells are initially 
populated with hydrogen and helium in a primordial ratio. 

When considering a uniform medium surrounding the 
source, the gas is given initial ionized fractions as deter- 
mined by the RecFast recombination code* (Seager et al. 
2000) for the appropriate cosmology and redshift. The ini- 
tial temperature of the gas in each shell is also determined 
by RecFast and each shell is initially set to be expanding 
with the Hubble fiow. 

We then proceed to evolve the thermal and ionization 
states of these shells forwards in time in a series of short 
time steps. During each time step we begin by computing 
the input spectrum of photons emitted by the central source 
(QSO, stars or both). Given this spectrum, we compute rates 
of ionization and heating in the innermost shell and solve for 
the evolution of its properties by integrating the appropri- 
ate set of differential equations as desribed below. The in- 
put spectrum is then attenuated by the optical depth of this 
first shell and used as input for the second shell. This pro- 
cess is repeated until the outermost shell is reached (which 
is chosen to be at sufficiently large radius that the radiation 
field is attenuated to close to zero at all times during our 
calculation). In addition to changes in temperature and ion- 
ization state, the density of each shell evolves as it expands 
or contracts due to any initial velocity and pressure forces. 
This approach is similar to those in other recent papers, e.g., 
Thomas & Zaroubi (2008). 

Our calculations of the ionization and thermal evolu- 
tion of each shell use the same input physics as the IGM 
evolution model of Benson & Bower (2010). The density of 
each ionization, n^, state in a given shell is then given by 



1^ 
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where for each atomic species H or He, i refers to their ion- 
ization state (i.e., i = 1 and 2 for H and H^, and i = 3, 
4 and 5 for He, He"*" and He'^"''), n^ is the number density, 
T is the temperature of the shell, V is the volume of the 
shell, Qi is the recombination rate for i (Verner & Ferland 
1996), Fe^i is the collisional ionization rate coefficient for i 
(Voronov 1997) and F^_i is the photo-ionization rate for i 
which is given by 



a',{E)ni 
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(2) 



where (t[ is an effective photo-ionization cross-section that 
accounts for the effects of secondary ionizations and is 



■* We use vl.4.2 of RecFast and include all of the modifications 
to the Hel recombination rate. 
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given by ShuU & van Steenberg (1985) (as re-expressed by 
Venkatesan et al. 2001): 



a'n{E) = (1 + 
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(4) 



where (^{E) is the actual cross section (Verner & Yakovlev 
1995) and 

0.4092n1.7592 



0Hi = 0.3908(1 - x" 
6He, = 0.0246(1 - X. 
teei = 0.0554(1 - X 



0.4049n1.6594 



0.4614n1.6660 



(5) 
(6) 
(7) 



In the above, S{E)dE is the number of photons emitted per 
second in the energy range E to E+dE by the central source 
and t{E; r) is the optical depth to radius r at energy E. 

Similarly, the evolution of the temperature of each shell 
is given by 






TdM ^ (s^-A^: 

^ di I fee "-tot 



(8) 



Here, 7 is the adiabatic index of the gas, S"^ is the rate 
of heating per unit volume due to all the heat sources (i.e. 
Compton heating and photo-heating) and A"^ is the rate 
of cooling per unit volume due to all the heat sinks (i.e. 
Bremsstrahlung cooling and various atomic processes), ritot 
is the total number density of atoms (H and He) and their 
ions per unit volume, T is the temperature of the shell and 
fee is Boltzmann's constant. 

In the above equation the first term represents adia- 
batic cooling due to the expansion of the shell. The second 
term accounts for the effects of changes in the mean atomic 
mass due to ionization and recombination processes. The 
final term accounts for the heating and cooling effects of 
the various processes that we now discuss below. 



Photoheating 

Photoionization heats the shell at a rate of 
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where Ei is the energy of the sampled photons which is asso- 
ciated with atom/ion number density rii, a' is the effective 
partial photo-ionization cross section (accounting for sec- 
ondary ionizations) for the ionization stages of H and He, 
n^(E) is the number density of photons of energy E, and Ei 
is the ionization potential of i. In the above, £ accounts for 
heating by secondary electrons and is given by (ShuU & van 
Steenberg 1985): 



S = 0.9971[1 - (1 - 



„0.2663\1. 31631 



(10) 



Compton Cooling/Heating 

Compton scattering of CMB photons from free electrons 
causes cooling or heating of the gas at a rate of (Peebles 
1968) 



Scompton ~ 4(7TaR (TcMb(1 + z)) 



4 Ticks 



where ctt is the Thompson cross section, aa is the radiation 
constant, Tcmb is the temperature of the CMB at z = 0, ric 
is the number density of electrons per unit volume and nic 
is the mass of an electron. 

For a typical source in our paper, we find that Compton 
heating is insignificant. The initial emission rate of ionizing 
photons for a lO'^ Mq starburst with a 10** M© BH (detailed 
in the next section) is ~ 1.3 x 10^^ photons s^^. The radius 
to which Compton heating is important (Ricotti et al. 2008) 
for this scenario at 2 = 10 is about 99 pc. As we will see, this 
is well below the 0.001-1 Mpc scales that are most relevant 
for I-front evolution and 21 cm signals in this work (§3); 
thus, Compton heating will not have a significant effect on 
our results. 
Single Electron Recombination Cooling 

Photon emission due to single electron recombination 
cools the shell at a rate 



Arec = T^bT [cVrH + 



(T)nH+ + arjj^+ (r)nHe 

+ "--Hc2+(^)"-Ho2+] Tie, 



(12) 



where a^ is the rate of the recombination processes for its 
respective atom/ion number densities, rii (Verner & Ferland 
1996). 

Dielectric Recombination Cooling 

Photon emission due to dielectric recombination cools 
the shell at a rate 

Adicicc = 40.74eV ad(T)nHe2+ne (13) 

where Od is the rate of the recombination process for He'^^ 
(Aldrovandi & Pequignot 1973; ShuU & van Steenberg 1982; 
Arnaud & Rothenflug 1985). 

Collisional Ionization Cooling 

CoUisional ionization leads to a cooling rate of 

Aion = [Eucei^{T)nu + EiicCi.i^„{T)nuc 

+ -£^Ho+"iH„+(^)"Ho+] "-0, 



(14) 



where Oi is the collisional ionization rate coefficient for the 
respective atom/ion of number density rii and Ei is the ion- 
ization potential of the respective atom/ion, H, He and He"^. 

Collisional Excitation Cooling 

Collisional excitation followed by radiative decay cools 
the shell at a rate: 



Aox = (ocoUh'^h + acoiijj^+ riHc+ ) ■ 



(15) 



where acoiiH and a;,oiiHo+ s-rs ths rates of collisional excita- 
tions involving H and He^ respectively (Scholz & Walters 
1991). 

Bremsstrahlung Cooling 

Finally, Bremsstrahlung emission cools the shell at a 
rate 



Ae 



16 / 2nkB \ 
3^3 [h^mij 



4c7reo 



c2Vr[7H+(T)nH+ 



+7He+ (r)riHc+ + 47ho2+ {T)n^i^2+] n^ 



(16) 



(rcMB(i + 2)-T),(ii; 



Here, eo is the permittivity of free space and 7 is the energy- 
averaged Gaunt factor (Sutherland 1998). 
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These coupled differential equations are solved numeri- 
cally using a standard Runge-Kutta method. 



3 RESULTS 

As noted earlier, we focus on early galaxies of typical mass 
~ 10^-10^" Mq in total mass and of approximate size a few 
kpc at most. We therefore perform most of our calculations 
at z — 10, with one calculation at 2 = 20 for comparison. 

To calculate the feedback from a typical 
QSO/starforming galaxy at these epochs, we compute 
the BH mass function at z = 10 using data that is pub- 
licly available from the Millennium Simulation database^ 
(Springel et al. 2005). In Figure 1, we show the computed 
BH mass function at z = 10, where we see that a typical 
quasar is powered by BHs in the mass range ~ 10^-10® Mq, 
which we use as a baseline for most of the cases considered 
in this paper. The turnover in Figure 1 may be partially 
due to the finite resolution of the simulation itself; in reality, 
we expect that the mass function should continue to slowly 
rise to somewhat smaller masses. In our models, the X-rays 
from the stellar populations are minimal, so we consider 
cases where the BH mass is typically 10^ M©, with some 
lower BH-mass cases (down to no BH) and one case with a 
BH mass of 10* Mq to derive an upper limit to the X-ray 
feedback. We assume that the duty cycle of the QSO is 100 
Myr for nearly all our cases but include one case with a 
low-mass BH QSO that has a shorter duty cycle of 10 Myr. 

Note that the the typical ratio of BH to stellar burst 
masses considered here are not consistent with the mea- 
sured ratio of the BH to stellar spheroid (bulge) mass of 
0.15% at z = (Giiltekin et al. 2009). Early galaxies differ 
from present day ones in that they must have a seed BH that 
grows with time over generations of starbursts and galaxy 
mergers. Today we measure the BH to star (or spheroid) 
mass ratio after these processes have happened but it is un- 
clear what this ratio would be for primordial galax;ies, or 
if this ratio remains constant down to lower galaxy masses 
(Greene et al. 2010). AGN observations indicate a possible 
lag in the peak of BH growth (and therefore AGN activity) 
relative to the peak in the star formation rate in early galax- 
ies, owing to gas dynamical effects between star formation 
and BH "feeding" (Hopkins 2011). There are additional un- 
certainties related to the gas fraction, the Eddington ratio 
etc. at high redshifts. Thus, wc provide a few example cases 
here but do not attempt to provide a cosmological sample 
of model galaxies. 

In order to distinguish the contributions of X-ray ion- 
ization relative to that from UV radiation, we consider three 
variations on each case with a starburst and QSO: one with 
the full spectrum including UV and X-ray photons from the 
source, one without the X-rays, and one with the X-rays 
alone. To do this, we need to define the boundary between 
what is considered an X-ray versus a hard UV photon, a 
quantity that has often not been clearly defined in the cos- 
mology literature on this topic (Chen & Miralda-Escude 
2008; Ricotti et al. 2005). At least some of this difference 



^ The Virgo-Millennium database is available at; http://www.g- 
vo.org/Millcnniuni/ 



arises from considering the spectrum at the source versus 
the emergent spectrum after processing through the gas in 
the galaxy. We choose 120 eV as the minimum threshold 
for what we consider an X-ray. This is consistent with the 
broader physics definition, but also with the impact of a 
typical X-ray on the IGM. We discuss this in detail in Sec- 
tion 3.3, but we note for now the well-known result that the 
mean free path (MFP) of X-rays varies substantially by X- 
ray energy. We show this explicitly in Figure 2: a 100 eV 
photon has a MFP of 0.1-0.2 Mpc whereas a 1 keV photon 
has a MFP that is larger by more than 3 orders of magni- 
tude. Note too the "ranking" of the three species in this plot 
- He I has the lowest MFP at all energies, representing the 
bottleneck for X-rays that results in secondary ionizations 
for H I (Section 2). 

3.1 Feedback from First Stars and QSOs 

We begin by examining a number of cases at z = 10 that in- 
volve varying combinations of starburst and BH masses. The 
plots all show cases with and without X-rays, and one with 
X-rays only (i.e. no lower energy photons). We begin with a 
lO'^ Mq starburst with a 10*^ Mq BH, hereafter referred to 
as the standard case. Figure 3 displays the ionization and 
temperature profiles as a function of distance from the cen- 
tral starburst/QSO source at z = 10, for the species H II, 
He II, and He HI. The red and green curves respectively 
show the evolution of the ionization and temperature curves 
at times 10 Myr and 100 Myr after the source turns on. The 
X-rays contribute from ~ a few percent up to full ionization 
in different H/He species at IGM scales (10-100 kpc), and 
heating of the order 10*-10^ K. Although the panels with 
and without X-rays (the upper two panels) look very similar 
at first glance, we note the extended tail of low-level ioniza- 
tion in H II and He II (but not He III) beyond the I-front: 
the signature of X-ray ionization. This can be seen in the 
red curves (10 Myr) on physical scales of 0.1-0.2 Mpc. 

We also consider cases where the BH mass and QSO 
duty cycle are varied. This reveals the various contributions 
more clearly, particularly that from X-rays. The results are 
shown in Figures 4, 5 and 6, where we can see that in- 
creasing (or decreasing) the BH mass or the duty cycle sim- 
ply "dials up" (or "dials down") the effects of ionization. 
For the higher BH mass, the X-ray I-fronts advance further 
and reach higher values of ionization. Nevertheless, the high 
temperatures of lO'' K and strong ionization effects from 
X-rays at large scales found by some authors, e.g. Thomas 
& Zaroubi (2008), are not reproduced here, possibly arising 
from differences in model assumptions and input spectra 
(discussed further in Section 3.3). 

Comparing the curves for the X-rays-only case for QSO 
BH masses of 0, 10^ Mq and 10* Mq, we see that X-rays can 
make a difference. Perhaps X-rays can become competitive 
with UV ionization only when the BH masses approach 10* 
Mq. Note that such high QSO BH masses are very rare at 
z = 10 (Figure 1), and likely nonexistent at z = 20 when 
the universe is younger and there has been little time to gain 
mass for a seed BH accreting at rates close to the Eddington 
value. Such 10* Mq or higher-mass AGN therefore may not 
contribute significantly to a cosmic X-ray background at z > 
10. Also, we point out that in all the figures the X-ray related 
features noted earlier (the tail of low-level ionization in H II 
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and He 11, but not He HI, at large radii) are evident in the 
upper two panels in each case. The exception is the case 
with only stars (10® M© starburst, BH mass of 0) where the 
figures with and without X-rays are (unsurprisingly) near- 
identical. 

Additionally, we ran cases with smaller masses in stars 
and BHs. One such case is shown in Figure 6, where the ion- 
ization and temperature profiles are displayed for a 10'' M0 
starburst with lO" Mq BH at z ~ 10, at times 10 Myr and 
100 Myr after the quasar turns on. Unlike previous figures in 
the paper, the no-Xrays case is not shown here, as it is very 
similar to the full spectrum case. The various panels show 
the curves for the full QSO spectrum (including UV/X-ray 
photons) and with X-rays only, with varying QSO duty cy- 
cles of 10 Myr, and 100 Myr. The ionization and maximum 
temperatures are lower over 10-100 kpc compared with our 
standard case but the role of X-rays for He I ionization is 
more clearly seen here than in most our cases, particularly 
in the X-rays only panel for a QSO duty cycle of 100 Myr. 

Other trends include variations with time or between 
species. Allowing the QSO/starburst source to be "on" for 
100 Myr advances the I-fronts for all cases and species rel- 
ative to the curves for 10 Myr, as expected. The tempera- 
tures, however, increase noticeably at 100 Myr only for the 
pure X-rays case; for the cases involving the full spectrum 
or without X-rays, the temperatures appear to saturate at 
a few tens of thousands of degrees Kelvin, and having the 
source on for longer timescales makes little difference. In 
addition, the He HI I-front mostly lags the H I-front but 
in some cases the He HI front almost catches up to the H 
I-front. Thus, it appears that these species' I-fronts can be 
coincident for sufficiently hard radiation. 

The He II ionization fraction exceeds that of H I by 
a small margin, particularly beyond the edge of the UV I- 
front. We recognize this as the characteristic tail of added 
secondary ionizations from X-rays, which manifest more 
strongly at larger physical scales where the UV photons do 
not penetrate as far. This can be seen best by comparing the 
no- X-rays and X-rays-only panels of all the figures in this 
section, where the He II front lags or is similar to the H I 
front when X-rays are absent but leads the H I front when 
only X-rays are present. This interplay between X-ray sec- 
ondary ionization and the ionization balance of H and He in 
the presence of hard radiation leads to ionization boundaries 
that are less sharp than in the UV-ionization case alone (see 
also Furlanetto & Oh 2008b on this point in relation to the 
morphology of helium reionization at lower redshifts, z ~ 
3). Last, in the case with only a 10® M© starburst (Figure 
5), we see that there is little difference between these two 
panels, as this case has low X-ray production. 

To test the variation with redshift, we perform the same 
calculations for our standard case assumptions at z — 20. 
Exploring redshifts lower than z ~ 10 marks the era of over- 
lapping I-fronts as reionization draws to an end, which our 
current treatment cannot model well. Additionally, there is 
not much H I remaining outside of galaxy halos to gener- 
ate an interesting 21 cm signal at the end of reionization, 
whereas the 21 cm signal is expected to be significant at 
z — 10-20. The calculations at 2: = 20 for our standard 
case are displayed in Figure 7 with the same three panels 
as in the ionization and temperature figures. As the IGM 
is denser and the recombination timescales are shorter, we 



show curves for times at 1 Myr and 10 Myr (rather than 10 
Myr and 100 Myr) after the source turns on. We see that 
the ionization curves at 10 Myr between the z — 20 case and 
our standard case at 2: = 10 have very similar shapes, with 
the 2 = 20 curves lagging the 2 = 10 curves, expected from 
the higher IGM densities at earlier times. Note however that 
the peak temperatures achieved in all of these cases remain 
similar, around 10"" K. 

We perform a simple estimate of the tradeoff between 
the local X-ray flux from a single galaxy versus the X-rays 
from a number of distant sources. The comoving number 
density of halos in our work with masses > 10* Mq is, n = 
1.147 (6.443 X lO""*) Mpc"^ at 2 =10 (20). This translates to 
an average spacing between such halos of ~ 0.95 (11.5) Mpc 
at 2 = 10 (20) . The emission rate of H-ionizing photons for 
a lO"' Mq starburst with a 10® M© BH (our typical case)® 
is S ~ 1.3 X 10''^ photons s^^. The associated X-ray photon 
production rate is ~ 1.3 x 10*^ (2.1 x 10'**) photons s"^ at 
300 cV and 1 keV respectively. If we assume a uniform IGM 
with no attenuation and that the visibility sphere for sources 
can go out to a maximum radius given by the MFP derived 
for X-rays as a function of energy in Fig. 2, then the critical 
distance from an individual galaxy source at which the flux 
of the source become equal to the background flux from 
sources of similar individual fluxes is 0.1-0.5 Mpc at 2 = 10 
for 300 eV to 1 keV X-rays. Thus, our results at 2 = 10, e.g. 
in Figure 3, could have additional contributions to X-ray 
ionization from neighboring galaxy halos at radii 0.1-1 Mpc, 
although this will be less of an issue at 2 = 20. In reality, 
we need to factor in realistic density profiles for the galaxies 
and the IGM, as well as the time variability of individual 
sources. We will pursue this in future work involving a full 
cosmological calculation through extensions to the current 
Galacticus code (see §4). 

Last, we note the oscillations in the He II fraction and 
temperature profiles in some of our models. We performed 
a number of checks to make sure these were not mere nu- 
merical effects. We found that these oscillations are robust 
to increases in the time resolution, ODE solver accuracy 
and number of radial shells used in our code. These os- 
cillations are also well-resolved radially, and have a near- 
constant wavelength, despite the logarithmically-spaced grid 
spacing in radius. What may be occurring is similar to the 
physics of the instability strip in stellar atmospheres. Inside 
the ionized region, the optical depth is very small, so the 
incident flux drops as 1/r^. The small H I, He I and He II 
fractions are determined by the balance between photoion- 
ization, coUisional ionization and recombination rates, while 
the temperature is controlled by the balance of photoheating 
and cooling rates. As we move outward in radius, this leads 
to a complex interplay between the photoheating rate, tem- 
perature and the He II fraction in the region of the He HI 
to He II transition, leading to the temperature and He II 
fraction oscillating with radius. This arises from our solving 
the time-dependent ionization and heating equations rather 



^ For comparison, S ~ 0.43 x 10^^ photons s ^ for a single 200 
Mq star in Chen & Miralda-Escude (2008), 5 ~ IQ-'''^ photons 
s~^ in Ricotti et al. (2005) (from the discussion related to their 
equation 4), and S ^ 10^''-10*'' photons s^^ for the BH mass 
range of lO'^-lO® Mq considered in Thomas & Zaroubi (2008). 
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than adopting the equilibrium solution. Given several of our 
idealized approximations here such as spherical symmetry, 
we do not expect this effect to have a significant impact, 
particularly on the 21 cm signal which we discuss next. 

3.2 Radio Signatures 

Over the last decade, there has been a growing literature on 
the 21 cm radio signals arising from the percolation of reion- 
ization, i.e., the growth of ionized bubbles around the first 
luminous sources and the associated heating (Zaldarriaga 
et al. 2004; Chen & Miralda-Escude 2004, 2008; Kuhlen et al. 
2006; McQuinn et al. 2006; Furlanetto et al. 2004; Furlan- 
etto & Pritchard 2006; Furlanetto et al. 2006; Pritchard & 
Furlanetto 2007; Thomas & Zaroubi 2008; Ripamonti et al. 
2008; Santos et al. 2008; Morales & Wyithe 2010). The sig- 
nature is expected to be absorption (emission) against the 
CMB if the ionized region is colder (warmer) than the CMB 
at those epochs. Forthcoming interferometric experiments at 
radio wavelengths, such as LOFAR and SKA, are predicted 
to be able to resolve ionized bubbles of size ~ 100 kpc up to 
a few Mpc. The dominant signal arises from the coupling of 
the spin temperature of neutral hydrogen with the kinetic 
temperature of the background IGM gas. After recombina- 
tion, the IGM cools as (1 -|- z)^ whereas the CMB cools as (1 
-I- z), leading to a 21 cm absorption signal from the neutral 
IGM gas. At later epochs, the spin states of hydrogen come 
into equilibrium with the CMB, leading to a decreasing 21 
cm signal. As the first stars and quasars turn on, a 21 cm 
emission signal is generated through coupling the spin states 
with the scattering of Lya photons and other processes. 

Here, we follow the formalism outlined in Chen & 
Miralda-Escude (2008). As we do not follow the detailed 
cosmological evolution of a distribution of ionized bubbles, 
we model the spin temperature of H I at a fixed redshift as: 



Ts 



TcMB + (j/a +yc)Tk 



(17) 



1 + J/« + J/c 

where Tcmb is the CMB temperature at that redshift {z — 
10 in our cases unless otherwise specified) and Tk is the gas 
kinetic temperature (which is a function of distance from 
the source). The y-coefficients are related to the coupling 
arising from Lya photons (j/q.) and from collisions (j/c). The 
coefficient t/c is taken from Chen & Miralda-Escude (2008) 
and Kuhlen et al. (2006). The coefficient ya is the Lya cou- 
pling term arising from the Wouthysen-Field effect. We use 
the expressions for y^, from Chen & Miralda-Escude (2008), 
Zaldarriaga et al. (2004), and Pritchard & Furlanetto (2007), 
with additional parameters from Hirata (2006). In the cases 
considered here, Lya coupling dominates over other terms 
such as coUisional coupling. We specifically include the Lya 
photons from the stars and/or QSO emission in our mod- 
els, as well as the auxiliary Lya photons arising from X-ray 
ionization (Chen & Miralda-Escude 2008; Vcnkatesan et al. 
2001; Shull & van Steenberg 1985). 

This leads to a brightness temperature (measured as a 
differential from the background CMB temperature at that 
epoch) given by: 



5Tb = 40 mK 




(18) 



When this calculated brightness temperature, STi,, lies above 
the CMB temperature at that epoch, the ionized region will 
be seen in emission against the CMB. Conversely, regions 
beyond the I-front that lie below the CMB temperature will 
be seen in absorption against the CMB. 

In Figures 8- 13, we show the temperature profiles with 
radius for the spin temperature and gas kinetic temperature 
relative to the CMB temperature which is constant at a fixed 
redshift. We also show the 21 cm brightness temperature 
profile and include the full spectrum case (X-rays and UV 
photons) and X-rays-only cases for each set of curves. These 
scenarios span most of the cases discussed in Section 3.1 
involving a combination of starburst and QSO/BH masses 
(most of which are at 2; = 10, with two cases at « = 20). 

Some broad conclusions that are common to all the 
cases whose 21 cm signatures are shown are as follows. First, 
the curves for the spin temperature are characteristically 
peaked around the location of the stalled I-front. The tran- 
sition from fully ionized within (with zero STi,) to the neu- 
tral IGM gas occurs beyond the I-front in each case, with 
peak values for Ts reaching ~ 10*-10^ K in our cases, and 
peak values for the STt emission signal around 30-40 mK. 
Negative 5Tb values, corresponding to an absorption signal 
relative to the CMB, occur on scales between 0.1 and 1 Mpc 
at 2: = 10 in our models and have low net values of ~ to 
a few mK, and larger values of ~ 20-60 mK on scales of 
0.01-0.1 Mpc at z = 20. We discuss this further below. 

Second, the curves in each case corresponding to the X- 
rays only case for each starburst/BH scenario consistently 
lag the curves for the corresponding full spectrum case. This 
is most dramatically seen in the stars-only case (Figure 13), 
a 10^ solar-mass starburst with no QSO/BH), where the 
X-ray production is low. Here, the maximum values of 5Tb 
occur between 1 and 10 kpc for X-rays only and at about 
50 kpc for the full spectrum. This case also reveals the in- 
herently "fuzzy" ionization fronts associated with X-rays, 
relative to the sharp I-fronts of UV radiation - note the 
gradual transition in spin temperatures for the X-rays-only 
case spanning nearly two orders of magnitude in scale. In 
contrast, the case of the 10^ solar-mass starburst with 10* 
solar-mass QSO/BH (Figure 12) reveals that the cases with 
and without X-rays barely differ in the location and peak 
values of Ts and 5Tb (emission in the latter). This arises 
directly in the strong contribution of X-rays to the overall 
ionization budget in this scenario. Ironically, it seems that 
the greater the X-ray production of a source, the less likely 
it is have a distinguishing X-ray-related signature at 21 cm. 

These results reveal one of the key goals of this paper: 
the difference in the topology of reionization between X-ray 
and UV ionization scenarios, and their impact on 21 cm 
predictions. Although X-rays do penetrate deeper into the 
IGM than do UV photons (leading to the moderate gains 
in ionization and temperature mentioned earlier), their "I- 
fronf's trail the UV I-fronts and therefore the UV-associated 
21 cm signal. This could therefore "blur" the signatures of 
the growth of ionized bubbles around first-light sources, and 
alter predictions for observing the percolation of reioniza- 
tion (see the semi- numerical simulations of Warszawski et al. 
2009 on this point). We note that a cosmological scenario in 
which X-rays alone are generated is not well-motivated phys- 
ically. Rather, the figures in this section show that the differ- 
ing scale-dependent ionization from X-rays and UV photons 
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lead directly to 21 cm signals that can be distinguished from 
each other. 

We consider time evolution in two cases for the same 
source of a 10^ solar-mass starburst with a 10® solar-mass 
QSO/BH: at z = 10 for times of 1 Myr and 10 Myr after 
the burst/QSO turn on (Figure 8 and Figure 9), and at 
z = 20 for times of 0.1 Myr and 1 Myr (Figure 10 and 
Figure 11). The main effects of the advancing I-front with 
time on the 21 cm signal at a fixed redshift are the following: 
a similar advancing of the spin temperature curve's peak, 
and therefore that of STb, from a few tens of kpc to about 
100 kpc, and, a decreased peak value in Ts- This is mostly 
due to the rapid falloff in the Lya flux at increasing radii 
(going as r~^), which leads to a decreased coupling between 
the gas H I and the source radiation. The important role 
of this LyQ photon coupling is manifested also through the 
slight increase in the positive values (emission signal) of STb 
and the increased negative values of STt (absorption signal 
at 21 cm) at 2; = 20 relative to z = 10, arising from the 
closer location of the I-fronts to the source with increasing 
redshift. These effects are discussed in more detail in the 
next section. 

3.3 Comparison with Other Works 

Here, we compare our results and model assumptions with 
those from papers in the recent literature addressing X-ray 
and/or helium ionization, and the resulting 21 cm signals. 
We find that our results are, for the most part, in agree- 
ment with the findings of other groups when we make similar 
model assumptions. We also comment on the theoretical as- 
sumption of passive X-ray production tied to star formation 
at high redshifts. 

In Chen & Miralda-Escude (2004) and Chen & Miralda- 
Escude (2008), the emergent spectrum is based on a radia- 
tive transfer calculation starting with a stellar blackbody 
spectrum. There is no stated definition to distinguish be- 
tween X-rays and UV radiation, so that (as in some works on 
this topic) it is unclear where the X-ray/UV photon bound- 
ary lies. To compare their results with ours, we started with 
the blackbody spectrum from equation 8 in Chen & Miralda- 
Escude (2008) . The range of Pop III star masses that they 
consider (25-800 Mq) leads to a relatively narrow range of 
blackbody temperatures, T^s ~ 1.06 x 10^^-1.17 x lO'^ K. 
In Figure 14, we show the blackbody energy output (the 
Planck energy density, in units of power per unit area per 
unit solid angle per unit frequency) as a function of energy 
for a 25 Mq and a 1000 Mq star. There is little difference 
between the two cases - essentially nearly all Pop III stars 
have the same energy output (Bromm et al. 2001; Tumlinson 
et al. 2003). 

However, what is relevant here for us is the cutoff be- 
tween UV and X-ray photons. The strict definition of X-rays 
has a lower limit of 120 eV for X-ray energies. In Figure 14, 
we see that the energy curves are relatively flat for energies 
of 20-40 eV and start to decline steeply above 100 eV. We 
were able to reproduce Figure 1 in Chen &: Miralda-Escude 
(2008) only for the X-ray threshold energy lying at about 
30 eV. Such "X-rays" can make a substantial addition to 
the UV-only ionization case, owing to the large numbers of 
photons below 100 eV. However, placing the cutoff at 100 
eV or higher (where the spectrum is down by a factor of ~ 



100 relative to the peak), leads to the results in our earlier 
ionization figures, where X-rays can have a significant (but 
not dramatic) impact on IGM ionization and temperature. 

Also, Chen & Miralda-Escude (2008) consider a single 
200 Mq star embedded in a galaxy, versus our treatment 
of a starburst and/or QSO as a point source in the IGM. 
The I-fronts in their work are therefore a factor of 10-20 
closer to the source than in our results, leading directly to 
a lower Lya flux in comparison at large scales in our work. 
Consequently, the 21 cm absorption signal induced by the 
LyQ photons in our calculations is weaker relative to that in 
Chen & Miralda-Escude (2008)'^ or, e.g., Thomas & Zaroubi 
(2008). This reduced signal is seen as a minor dip, rather 
than a larger trough, in the 21 cm brightness temperature 
beyond the I-front location in the right panels of the fig- 
ures for the z — 10 cases in Section 3.2. We have checked 
this by artificially placing the I-fronts in our cases in the 
21 cm calculations at closer radii (by ~ a factor of 10) and 
are able to reproduce the 21 cm brightness temperature ab- 
sorption signal of Chen & Miralda-Escude (2008) and other 
works. Note that for the z = 20, t = 0.1 Myr case (Figure 
10), the absorption trough becomes more noticeable as the 
I-front has not advanced as far. This verifies the critical role 
of the invere square dropoff of the Lya photon flux with dis- 
tance from the source for 21 cm absorption (discussed in an 
earlier section). Note also that our results for the predicted 
amplitudes of the spin temperatures and the 21 cm emission 
signal are in agreement with other papers in the literature. 

In Bolton et al. (2009), McQuinn et al. (2009), Furlan- 
etto & Oh (2008b) and Furlanetto & Oh (2008a), the authors 
focus on helium reionization by quasars at 2 ~ 3. Furlanetto 
& Oh (2008a) do not include the effects of X-ray heating in 
their calculations of helium ionization. Bolton et al. (2009) 
find a relatively modest gain in IGM temperature (of order 
lO** K) resulting from hard radiation, partly owing to the 
heating in underdense parts of the IGM (particularly fossil 
He HI regions) achieving their maximal heating early on in 
the process of reionization (see also Venkatesan et al. 2003 
on this point). This maximum IGM temperature of ~ 10'* 
K (comparable to the results of McQuinn et al. 2009) lies 
within the range of our findings, with the caveat that at 
z ~ 3 the IGM is far less dense than at z ~ 10-20, and, 
additionally, the IGM hydrogen is completely reionized at 
z ~ 3, freeing up some of the UV photons and secondary 
electrons from He I ionization. We also compared our results 
with Kuhlen et al. (2006) - these authors do not include Lya 
coupling in their 21 cm calculations but we approximately 
reproduce their results on spin temperature values. 

Thomas & Zaroubi (2008) have examined the feedback 
from early stellar populations and quasars, and the asso- 
ciated 21 cm signature. We found that we were unable to 
reproduce many of their results, including the high level of 
X-ray heating (T > lO^-lO" K) as well as the results in 
their Figures 12-13. Some of this may arise from incomplete 
models of the high-z galaxy distribution and that the stellar 
spectra have been simplified as a blackbody source. We do 



^ See these authors' discussion in Sec. 2 of tiicir paper of tiic typ- 
ical size of Lya spheres in their work being a few tens of kpc, and 
of their assumption that the fraction of X-ray energy converted 
to Lya photons is 100%. 
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however find somewhat similar IGM temperatures and 21 
cm brightness temperature values close to those computed 
in Ripamonti et al. (2008), although these authors focus on 
X-rays from BHs in the pre-reionization IGM. A direct com- 
parison is challenging as we do not have a full cosmological 
calculation in this paper with a halo distribution function 
characterizing the ionization feedback. 

Last, Ricotti et al. (2005) consider the formation of a 
strong X-ray background generated at z > 10 by early black 
holes, with the specific aim of explaining the WMAP results 
of a (at the time high) electron scattering optical depth. In 
their work, the X-rays are created by accretion onto "seed" 
black holes that are assumed to have formed from earlier 
generations of Pop III stars. Hence, the total X-ray emis- 
sivity is proportional to the total mass in such black holes 
in high-z halos, which is in turn proportional to the total 
mass in Pop III stars. That is, the production rate of X-rays 
is tied effectively to the star formation efficiency in high- 2; 
galaxies (through the black hole accretion rate). This is an 
assumption made in a number of papers, e.g., Santos et al. 
(2008) who examined the role of inhomogeneous X-ray and 
LyQ radiation fields at 2 >10 for 21 cm signatures of H reion- 
ization. However, it remains to be seen how well this series 
of connections hold at the low black hole masses anticipated 
in the first galaxies (e.g., do these low-mass BHs even ac- 
crete at the Eddington rate?). The Magorrian relation may 
not hold at low to moderate BH masses in galaxies (Greene 
et al. 2010), making the scaling of X-ray production with 
star formation rates and BH masses more ambiguous at low 
BH masses. 

Ricotti et al. (2005), like Chen & Miralda-Escude 
(2008), do not explicitly distinguish between soft X-rays and 
hard UV photons in their calculations. In these works, the 
boundary between UV and X-ray photons is related to the 
local column density of absorbers and the emergent power 
spectrum after radiative transfer, with the column density 
being a free parameter. Another important related param- 
eter is the escape fraction of ionizing radiation, /esc, for X- 
ray and UV photons, which is effectively calculated locally 
through the emerging fiux at each radius (or cell) in our work 
and for the above papers. Variations in parameters such as 
the local absorber column density, /osc and reduced gas den- 
sities within galaxies owing to feedback effects could harden 
the source spectrum within the ionization bubbles, leading 
to potentially higher temperatures than we have found here. 
Although there is no straightforward way within the scope 
of our semi-analytic work to directly reproduce local /osc 
and density-feedback effects from numerical simulations, we 
mention these caveats and note that we are able to reproduce 
the results of Chen & Miralda-Escude (2008) by lowering the 
X-ray/UV boundary to 30 eV or placing the I- fronts closer to 
the source, both of which effectively harden the local ionizing 
spectra. Last, Ricotti et al. (2005), Chen & Miralda-Escude 
(2008), and Thomas & Zaroubi (2008) have a fully cosmo- 
logical calculation that keeps track of the evolving spectra 
of stellar and QSO populations. Therefore, soft X-rays from 
high- 2; sources are redshifted and can become important for 
hard UV ionization at later epochs. Our current results do 
not factor in this effect, but we plan to extend this work 
in the near future to fully cosmological calculations that in- 
clude realistic galaxy profiles and the redshift evolution of 
galaxy halos and their radiation fields. 



4 CONCLUSIONS 

We have examined the effects of X-rays from high-redshift 
quasars and stars when acting in combination with hard 
UV ionizing radiation from these sources. We find that, rel- 
ative to hard UV radiation. X-rays may not dominate the 
ionization and thermal history of the IGM, and contribute 
modest increases to the IGM ionization at 2 ~ 10 and con- 
tribute of order 10^-lO^K to the IGM temperatures. This is 
in contrast with some earlier works in which X-rays could 
cause IGM heating up to 10^ K and near-total reionization 
at 2=10-20. While some of this may be due to our simpli- 
fied models, we believe that most of the difference between 
our results (where we include the X-rays coming from indi- 
vidual sources) , and those of other works deriving high IGM 
temperatures and ionization from X-rays at z = 10-20, arise 
from the latter's assumption of strong X-ray production that 
is tied to the star formation rate at high redshift. 

We also examined the 21 cm signatures of various cases 
involving combinations of stars and black hole masses, and 
find that the 21 cm signal of X-ray versus UV ionization 
could be distinct, resulting from their differing contributions 
to the topology of reionization. We find that the brightness 
temperature emission expected from X-rays alone occur at 
smaller scales than that from UV radiation. The different 
spatial scales at which they manifest may therefore "blur" 
the 21 cm signature of the percolation of reionization around 
early halos, depending on whether a cosmic X-ray or UV 
background built up first. An X-ray background may not 
significantly precede a UV background, as a typical X-ray 
photoionization timescale exceeds the Hubble time for 2 > 
10. From our simpified treatment, it is unclear whether there 
is a cosmological epoch when the IGM's thermodynamic and 
ionization properties are determined mostly by X-rays. The 
role of X-rays versus hard UV radiation can also be tested 
through their interactions with the CMB, where the rela- 
tive strengths of their contributions to reionization as well 
as the redshifts that they dominantly contribute at can be 
constrained through the CMB polarization power spectrum 
at large angular scales. The currently operating all-sky CMB 
mission Planck may be able to distinguish such scenarios. 

For sufficiently hard radiation from sources, the H II 
and He HI I-fronts may lie very close to each other. Although 
our calculation is ID in nature, this result will impact the 
escape fraction of ionizing radiation from primordial galax- 
ies, and the geometry of bubbles and chimneys as ionization 
proceeds from these galaxies. We hope to examine this prob- 
lem in a future work. 

To further explore the evolution of first-light sources 
and the IGM with redshift, we will extend our current cal- 
culations using the Galacticus code to a fully cosmologi- 
cal framework that includes evolving dark matter halos, 
galaxy/BH formation, and evolving stellar/QSO popula- 
tions with time-dependent radiation fields. This will permit 
a self-consistent calculation of IGM reionization, and allow 
us to derive predictions for the growth and evolution of a 
cosmologically representative distribution of ionized bubbles 
as a function of redshift. We can also calculate the bub- 
bles' thermal properties, as well as the statistical proper- 
ties of the bubble population, such as the mean size of ion- 
ized and neutral regions and power spectra of 21 cm emis- 
sion or absorption relative to the CMB (utilizing the known 
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correlation properties of the dark matter halos which host 
the sources). Such predictions will be tested by data from 
CMB space telescopes such as Planck, and ground-based ra- 
dio telescopes that are designed to map the percolation of 
reionization around first-light sources. These observations, 
coupled with our detailed theoretical predictions, will addi- 
tionally place strong constraints on the populations of ioniz- 
ing sources at intermediate to high redshifts and, therefore, 
on the properties of early generations of galaxies and AGN. 
The resulting improvements in our understanding of these 
early objects will permit more robust predictions to be made 
for other observing programs, such as those of the James 
Webb Space Telescope, which will probe similar galaxy/QSO 
populations. 
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Figure 1. The BH mass function at 2 = 10, using public data 
from the Millennium Simulation database. Note the peak around 
10^—10® M0. See text for more discussion. 
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Figure 3. The ionization and temperature profiles for a 10^ solar- 
mass starburst with 10® solar-mass QSO/BH at 2 ~ 10 (our stan- 
dard case). The solid, dotted and dashed lines represent the frac- 
tion of H II, He II and He III respectively. Red and green curves 
show the curves at times 10 Myr and 100 Myr after the quasar 
turns on. The upper left and upper right panels display the cases 
with the full QSO spectrum with UV photons, but that exclude 
and include X-rays from the central QSO. The lower panel shows 
the effects arising from X-rays alone. 



Figure 2. The mean free path in Mpc for H I, He I, and He II 
at 2 = 10 for photon energies ranging from 0.1 to 1 kcV. 
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Figure 4. The ionization and temperature profiles for a 10^ solar- 
mass starburst with 10* solar-mass QSO/BH at 2 ~ 10. 



Figure 5. The ionization and temperature profiles for a 10® solar- 
mass starburst only (no QSO/BH) at z ~ 10. 
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Figure 6. The ionization and temperature profiles for a 10^ solar- 
mass starburst with 10* solar-mass QSO/BH at 2 ~ 10. The left 
and right panels in each row display the cases with the full QSO 
spectrum (including UV/X-ray photons) and with X-rays only. 
The upper row has a QSO duty cycle of 10 Myr, and the lower 
row has a QSO duty cycle of 100 Myr. Red and green curves show 
the curves at times 10 Myr and 100 Myr after the quasar turns 
on. Unlike previous figures in the paper, the no-Xrays case is not 
shown here, as it is very similar to the full spectrum case. 



Figure 7. The ionization and temperature profiles for a 10^ solar- 
mass starburst with 10^ solar-mass QSO/BH at 2 = 20. The 
legend is the same as in earlier figures. See text for explanation. 
Curves are shown for timescales of 1 and 10 Myr (rather than 
10 and 100 Myr as in all preceding figures), owing to the shorter 
IGM recombination timescales at 2 = 20 relative to 2 = 10. 







Figure 8. Left panel shows the temperature profiles with radius 
for the spin temperature (black curves), kinetic temperature (blue 
curves) and the CMB temperature (purple line). Right panel dis- 
plays the 21 cm brightness temperature profile. All cases are for a 
10* solar-mass starburst with 10^ solar-mass QSO/BH at 2 = 10, 
at times of 1 Myr after the burst/QSO turn on. In each case, the 
solid lines are for the full spectrum case (X-rays and UV radia- 
tion) and the dashed lines are for X-rays only. 
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Figure 9. The temperature profiles (left panel) and 21 cm bright- 
ness temperature (right panel) are shown for a 10"" solar-mass 
starburst with 10^ solar-mass QSO/BH at z = 10, at times of 
10 Myr (rather than 1 Myr) after the burst/QSO turn on. The 
legend is the same as for Figure 8. 




Figure 10. The temperature profiles (left panel) and 21 cm 
brightness temperature (right panel) are shown for a 10^ solar- 
mass starburst with 10® solar-mass QSO/BH at z = 20, at times 
of 0.1 Myr after the burst/QSO turn on. The legend is the same 
as for Figure 8. Note the deeper trough in 5T^, (relative to the 
same case at z = 10) in the absorption signal against the CMB 
at scales of ~ a few tens of kpc. 
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Figure 13. The temperature profiles (left panel) and 21 cm 
brightness temperature (right panel) are shown for a 10® solar- 
mass starburst with no QSO/BH at z = 10, at times of 1 Myr 
after the burst turns on. The legend is the same as for Figure 8. 
The difference is more pronounced between the cases with only 
X-rays versus X-rays and UV radiation, owing to the low X-ray 
production of stars in our models. 




0.1 1 10 100 1000 

Energy (eV) 

Figure 14. A comparison of the blackbody energy output (the 
Planck energy density, in units of power per unit area per unit 
solid angle per unit frequency) from a 25 Mq and 1000 A/q star. 
Note the relative flatness of the curves at energies of 20—40 eV; 
beyond 100 eV the curves decline steeply. 



Figure 11. The temperature profiles (left panel) and 21 cm 
brightness temperature (right panel) are shown for a 10^ solar- 
mass starburst with 10® solar-mass QSO/BH at z = 20, at times 
of 1 Myr (rather than 0.1 Myr) after the burst/QSO turn on. The 
legend is the same as for Figure 8. 





Figure 12. The temperature profiles (left panel) and 21 cm 
brightness temperature (right panel) are shown for a 10^ solar- 
mass starburst with 10* solar-mass QSO/BH at z = 10, at times 
of 1 Myr after the burst/QSO turn on. The legend is the same as 
for Figure 8. 



